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The microscopic spectral density of the QCD Dirac operator at nonzero baryon chemical potential 
for an arbitrary number of quark flavors was derived recently from a random matrix model with the 
global symmetries of QCD. In this paper we show that these results and extensions thereof can be 
obtained from the replica limit of a Toda lattice equation. This naturally leads to a factorized form 
into bosonic and fermionic QCD-like partition functions. In the microscopic limit these partition 
functions are given by the static limit of a chiral Lagrangian that follows from the symmetry breaking 
pattern. In particular, we elucidate the role of the singularity of the bosonic partition function in 
the orthogonal polynomials approach. A detailed discussion of the spectral density for one and two 
flavors is given. 
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I. INTRODUCTION 



One of the features that makes QCD at nonzero baryon chemical potential both elusive and interesting is that 
the Euclidean Dirac operator does not have any Hermiticity properties. The non-Hermiticity occurs in an essential 
way because a finite baryon number density is obtained by promoting the propagation of quarks in the forward time 
direction and inhibiting the propagation in the backward time direction. Because of the non-Hermiticity the quark 
determinant attains a complex phase, which prevents the analysis of the QCD partition function by means of proba- 
bilistic methods. An approximation that is often used at zero chemical potential is to ignore the fermion determinant 
altogether. However, this approximation fails dramatically at nonzero chemical potential El- In particular, the crit- 
ical chemical potential at zero temperature is determined by the pion mass instead of the nucleon mass. Therefore, 
the main source of information for QCD at nonzero baryon density is based on simplified models or on perturbative 
expansions at very high densities which may never be accessed experimentally. Lattice QCD simulations are feasable 
in a region where the ratio of the chemical potential and the temperature is sufficiently small so that extrapolation 
from an imaginary chemical potential [2j, |3| or from zero chemial potential |^ becomes reliable (for a critical review 
of the field and additional references, we refer to 0). 

In this paper, we will investigate an observable for which nonperturbative results can be obtained for unquenched 
QCD at nonzero baryon density. This observable is the spectral density of the QCD Dirac operator. At nonzero 
chemical potential its support is a two-dimensional domain. The naive argument that observables do not depend on 
the chemical potential at zero temperature does not apply to the spectral density of the Euclidean Dirac operator. 
This puzzling fact can be understood as follows. To define an eigenvalue density we need both the eigenvalues and 
the complex conjugate eigenvalues or the Dirac operator and its complex conjugate. However, complex conjugation 
is equivalent to changing the sign of the chemical potential. Therefore, the generating function for the QCD Dirac 
spectrum has to contain both quarks and conjugate quarks in the valence sector that have an opposite baryon charge. 
This opens the possibility that the low energy limit of this theory contains Goldstone bosons made out of quarks 
and conjugate anti-quarks which have a nonzero baryon number H 0- More formally, the chemical potential in 
the generating function breaks the flavor symmetry resulting in a chiral Lagrangian that becomes dependent on the 
chemical potential. The region of the Dirac spectrum we wish to analyze determines the mass of the valence quarks 
d In order to determine the eigenvalue density of the Dirac operator at small eigenvalues the masses of the 
valence quarks are equally small. Consequently, the mass of the charged Goldstone modes is small and the generating 
function depends on the chemical potential even at small values of fi. 

The low energy effective theory for the generating function of the Dirac spectrum is a theory of Goldstone modes 
determined by spontaneous breaking of chiral symmetry. This effective theory is a version of a chiral Lagrangian 
[To| which takes into account the presence of the conjugate quarks. The allowed terms in the effective theory are 
determined by the flavor symmetries and the way they are broken by the quark masses and the chemical potential: 
the effective theory must break the flavor symmetries in precisely the same way as in the generating function we 
started from. To be specific, the Goldstone field, U, takes values in the coset U £ SU(Nf + 2n), where n is the 
number of additional quark - conjugate quark pairs. The leading order chiral Lagrangian in the standard counting 
scheme is the usual non- linear sigma model [TT| . The dependence of the chemical potential is completely fixed by the 
flavor symmetries. It enters as the zeroth component of an external vector field through a shift of the Euclidean time 
derivative (see e.g. El El El El EI) 

d U ^d U-i4B,U], (1) 

where /i is the chemicalpotential and B is the charge matrix. Since the baryon charge of conjugate quarks is opposite 
to the standard quarks 6] the charge matrix B is not proportional to the identity and the commutator is non- vanishing. 

Although the eigenvalue spectrum of the Dirac operator is not a direct physical observable, we hope that detailed 
knowledge of this observable in a nonperturbative domain where no other information is available, will ultimately lead 
to a better understanding of the problems that hinder numerical simulations at nonzero baryon chemical potential. 
At zero temperature these problems are manifest unless fj?F 2 V <C 1, where F„ is the pion decay constant. In this 
paper, we will focus on a scaling regime where the product \j?F 2 V is fixed as the volume, V, is taken to infinity. The 
fixed combination /j, 2 F 2 V can take any value and our results will show the effect of the fermion determinant on the 
Dirac spectrum. The quark masses will be taken such that scales like fj, with m^F^V kept fixed. To be specific, 
we consider the range [Slllq 

Wtt, ,u<C-<A, (2) 



as the volume L 4 — V is taken to infinity, and A is the scale of the lightest non-Goldstone particle. This regime is 
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sometimes known as the ergodic domain or as the epsilon regime 1 of chiral perturbation theory p"?l ]. In this domain 
the zero modes of the Goldstone fields dominate the partition function which reduces to a static integral over the 
Goldstone manifold [T3 . ITtT ]. The static integral is completely determined by the flavor symmetries of the QCD 
partition function. This implies that any theory with the same flavor symmetries and flavor symmetry breaking is 
described by the same static integral. The only memory of the underlying theory is two coupling constants, namely 
the chiral condensate E and the pion decay constant F^. Hence the partition function in this limit is universal in 
the sense that all theories with a given symmetry breaking pattern and a mass gap will have a partition function 
given by the same static integral over the Goldstone manifold provided that m, , // <C 1/L < A [Pj, Il9l |20j . The 
simplest theories in this class are invariant random matrix theories in the limit of large matrices. Because of the large 
invariance group it is sometimes simpler to analyze the random matrix theory rather than to evaluate the integral 
over the Goldstone manifold directly. 

In order to derive the microscopic eigenvalue density of the Dirac operator, which is defined by rescaling its 
eigenvalues as ZkVY>, we start from the microscopic partition function with Nf + In flavors. Then we take derivatives 
with respect to the mass of the additional n quarks and n conjugate quarks, and finally we remove the additional 
flavors by taking the limit n — > 0. This procedure is known as the replica trick [2lJ. Since the microscopic partition 
functions are only known for integer n there is no guarantee that a correct nonperturbative answer is obtained this 
way. After a two decade long discussion two closely related methods have emerged that result in the correct 
nonperturbative results: the replica limit of the Painleve equation pit , |24| and the replica limit of the Toda lattice 
equation [25l |26( . Both the Painleve equation and the Toda lattice equation are well-known in the theory of exactly 
solvable systems (see for example [27], y8|). The Painleve equation is a complicated nonlinear differential equation, 
whereas the Toda lattice equation is a simple two step recursion relation. For that reason, it is much simpler to 
work with the Toda lattice equation. The advantage of the Painleve equation is that nonperturbative results can be 
obtained from fermionic partition functions only |23l . Recently, the consistency of the replica limit of the Toda lattice 
equation [29| and the supersymmetric method f30j was established. A necessary ingredient for the applicability of the 
Toda lattice method is that the QCD partition function satisfies the Toda lattice equation. This is the case if the 
QCD partition function is a r- function [2^|. For the low-energy limit of QCD at /x = this was shown for fermionic 
partition functions in [3ll 132] ] . These results were extended to bosonic and supersymmetric partition functions as well 
as to QCD at nonzero chemical potential in l26l |29| . 

In the replica limit of the Toda lattice equation, generating functions with bosonic quarks appear in addition 
to generating functions with fermionic quarks. While the fermionic generating functions for the spectral density at 
nonzero baryon chemical potential are known from the integral over the Goldstone manifold [26l l33l | . only the simplest 
case with bosonic quarks is known ptjj . The relevant supersymmetric generating functions, for which the Goldstone 
manifold is a supermanifold, have proved to be quite challenging to compute, and no explicit expressions have been 
obtained so far. In the present paper we derive these generating functions from the random matrix model introduced 
in [34|]. The flavor symmetries and their spontaneous and explicit breaking are exactly the same in the large N limit 
of this random matrix model and in full QCD. Since both theories have a mass gap, in the microscopic limit, they 
are both given by the same integral over the Goldstone manifold. 

In the framework of random matrix models, spectral correlation functions can also be computed by means of 
(bi- ) ort hog onal polynomials in the complex plane. The complex orthogonal polynomial approach was developed in 
[35L l36l l37l l38l \3^L |40| and was applied to unquenched [34[ QCD Dirac spectra at nonzero chemical potential. The 
outcome of our analysis here is that both the Toda lattice and the orthogonal polynomials approach produce equivalent 
results. The divergence of the bosonic generating functions also sheds some light on universality of non-Hermitian 
random matrix models in general. 

Spectra of the non-Hermitian Dirac operator at nonzero baryon density have been obtained from lattice QCD in 
the quenched case 4lll42ll43| and for QCD with two colors j4JjE3 an d have been compared successfully to random 
matrix theory 0, ^ . l43t l45| . The microscopic spectral density of quenched lattice QCD at nonzero baryon density 
was first analyzed in 43j where quantitative agreement with analytical predictions |26ll37| was found in an asymptotic 
domain where the results derived from the chiral Lagrangian |26j agree with the expression obtained in |37j . The two 
flavor phase quenched partition function which is calculated in this paper does not suffer from a sign problem either 
and could be simulated on a lattice (see 0, ^3 f° r recent results) . A first analytical prediction for the unquenched 
QCD Dirac spectrum at nonzero baryon chemical potential in the non-perturbative regime was obtained in [34j. These 
results and the work presented in this paper could in principle be compared to numerical simulations in the region 
where the sign problem sets in. In particular, any proposal for a solution of the sign problem on the lattice can be 



1 The epsilon regime is the regime where m„- ~ 0{e 2 ) and 1/L ~ O(e). These conditions are more strict than the inequality For 
example, we are still in the ergodic domain if 1/L ~ 0(e 3//2 ). 
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tested against these predictions. 

For the spectral density with one dynamical flavor we show both analytically and numerically that large fluctuations 
occur for increasing value of the chemical potential. These findings further illustrate the difficulties encountered when 
one wishes to reproduce these results in lattice QCD simulations. 

In section [H] we explain the derivation of the spectral density of the QCD Dirac operator within the replica 
framework. The required generating functions are all given in terms of group integrals that represent low-energy 
effective QCD-like partition functions. The bulk of this paper is devoted to the calculation of the group integrals 
corresponding to bosonic quarks using the large N limit of a random matrix model. After introducing the relevant 
random matrix model for QCD with chemical potential in section HTT1 we compute the necessary generating functions 
by means of complex orthogonal polynomials in section lTVI In section[3we show that these partition functions satisfy 
the Toda lattice equation. This allows us to calculate the spectral density from the replica limit of the Toda lattice 
equation. We show in general terms that the results agree with a direct computation using orthogonal polynomials 
|34| . In section IVT1 we present our results for one, two and any number of dynamical flavors and compare them to the 
quenched and phase quenched results. It also contains a discussion of the thermodynamic limit of our results. This 
section is self contained and can be read independently of the previous sections. Our main findings are summarized 
in the conclusions, and several technical details are worked out in two appendices. 



II. THE SPECTRAL DENSITY FROM GENERATING FUNCTIONS 



At nonzero baryon chemical potential the Dirac operator is non-Hermitian so that its spectrum has a two- 
dimensional support in the complex plane. In this section we start with the definition of the spectral density for 
complex eigenvalues. Then we show how it can be obtained from a generating function with additional flavors. 

We consider the eigenvalues of the Dirac operator given by the eigenvalue equation 

(D v j v + ^ )ipj = Zjipj, (3) 

Since (D^7 r) )^ = — -D,,7, ; and (n'joY — M7o the Dirac operator is neither Hermitian nor anti-Hermitian and the 
eigenvalues are complex. Because of the axial symmetry {Z? ?? 7 ?? + /170,7s} = the nonzero eigenvalues come in 
pairs with opposite sign, ±Zj. The density of eigenvalues in the presence of Nf flavors is defined as the vacuum 
expectation value over a sum of delta functions in the complex plane at the positions of the eigenvalues, vanishing 
and non- vanishing, 

p^(z,z*,{ mf }^) + ^ 2 \z)^/^2S 2 (z-z 3 )\ , (4) 

\ i I QCD, v 

where fi is the chemical potential, Nf is the number of flavors and {rrif} — mi, . . . ,mN f are the quark masses. For 
later convenience we have not included the zero eigenvalues in our definition of the spectral density. The average over 
gauge fields in a fixed topological sector, u, is defined by 

^ J[d4 O rfe det(A,7„ + M7o + m f ) e^™^ 
QCD ^ f[dA]„ Ujii det(A, 7 r, + M7o + m f ) e-^(A) 

Since the phase of the product of the fermion determinants is non- vanishing, such expectation values of operators will 
in general be complex. In particular we will find that the unquenched spectral density is complex rather than real 
non-negative. 

The two-dimensional (^-functions appearing in the definition of the eigenvalue density can be obtained from the 
derivative 

(y^S^z-zA = ld„G? f {z,z*,{m f y,ij), (6) 

\ 3 I QCD, u 

with the resolvent defined by the average 
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One way to calculate the resolvent is by means of the fermionic replica trick. In the usual replica trick the resolvent 
is calculated from the identity 



E 



z 

3 



— = lim -d z T\{z + Zj ) n . (8) 



The problem is that, in most cases, the average of the right hand side can only be obtained for nonzero integer values 
of n, and the limit n — > can only be taken after a proper analytical continuation in n. For a non-Hcrmitian operator 
this procedure does not work. Instead, a modified identity has been suggested [I^ 

E— ^= Hm-&ni(2 + z;)| 2n - ( 9 ) 

3 3 

At the perturbative level the correct resolvent can be derived this way |48| . However, for nonperturbative calcula- 
tions this relation only leads to correct results if it is used in combination with the Toda lattice equation [25| . 
The generating function for the resolvent is thus given by 



■Zf" n ({m,}, z, z*; fi)= f [dA] v f[ dct(D nJn + M7o + m f ) | det(D, 7 , + ^7o + z)\ 2n e - s ™^. (10) 
The resolvent is obtained by 

G„ f (z,z*,{rrif};fj,) = lim -<9 Z log Zu I,n ({mf}, z, z*; fj,), (11) 



n^O n 

and the spectral density is given by 



p^(z,z*,{m f };n) + v5W(z) = lini—d z ,d z logZ^' n ({m f },z,z*;fj,). (12) 

n— 7T7J 

This construction has been used previously for perturbative 0, an d nonperturbative [24l l2fil | studies of the eigen- 
value density in quenched theories (i.e. for Nf = 0). 

The interpretation of the partition function (|l(Jf> is 0, Q that in addition to the usual quarks we have n quarks 
with mass z and n conjugate quarks with mass z*. Since 

det(D r) 7r; + M7o + z)* = det(7 5 (-Z? t? 7^ + fij + z*)~f 5 ) = det(D J) 7 J) - ^70 + z*), (13) 

conjugate quarks carry the opposite baryon number of quarks. The baryon charge matrix in the partition function 
(IIOII is therefore not proportional to the identity and the flavor symmetry is broken by the chemical potential. The 
specific way in which the flavor symmetries are broken determines the low-energy effective theory. In particular, the 
effective theory will depend on \i for energies where the pions dominate the free energy. Physically this is clear because 
mesons made out of quarks and conjugate anti-quarks carry a nonzero baryon number. The general prescription for 
constructing a chiral Lagrangian is to impose the transformation properties of the microscopic theory on the low- 
energy effective theory |l0j. For an external vector field the microscopic partition function is invariant under local 
gauge transformations of this field. Therefore, the derivative in the chiral Lagrangian has to be replaced by the 
corresponding covariant derivative (if! il2^ . 

In this paper we calculate the eigenvalue density in the microscopic limit for any number of flavors Nf, where 
rrifVT,, zVE, and /j, 2 F%V are held fixed as V — * 00. In this limit the kinetic terms factorize from the partition 
function and the low energy limit of the generating function (|10fl is given by [l4| 

Z? f ' n ({m f },z,z*;n) = f dU det(U) v e -^ 2 ^W,B][u-\B] + ^VTrM { u+u^) (w) 

JUeU(N f +2n) 

The quark mass matrix is given by M=diag(rni, . . . , mjVf 5 {z}n, { z *}n) and the charge matrix B is a diagonal matrix 
with 1 appearing Nj + n times along the diagonal and —1 appearing n times. An explicit expression for this integral 
was derived in |2(1 |33j . In this paper we will show that the partition functions (|14fl satisfy the Toda lattice equation 

™ tt ° g " ({m f },z,z,^- 2 (zz) [Z ^ {{mf },z,z*;^ ' ( } 



5 



where S z — zd/dz and 5 Z * — z*d/dz* . Since the replica limit of the right hand side of this equation is the spectral 
density 2 we find the remarkably simple expression l26l |29| 



Jfff, i m \.,A Zv ({m f },z,z ;^)Z l/ t ({m f }\z,z ; /i) 

p u (z,z ,\mfi,n) - — r ryN f ,n— o / r ; ^ ■ 

[Z v * ({mf};n)]' 



In (|16l) and elsewhere in this paper it is our notation to separate the fermionic and bosonic quark masses in the 
argument of the partition function by a vertical line. In addition to the fermionic partition functions, we will need 
to evaluate Z v f ' n for n = — 1, i.e. for one bosonic quark and one bosonic conjugate quark. In the quenched case 
(Nf = 0), this partition function is given by the integral pfj Eflj 



ZT-Hz,Z*',li) = limCe / f!, 0(U) e -^F^[U,B][U-\B] + iV^TrM^iU-IU-^)^ 

J der M U 

where dU6(U)/det 2 U is the integration measure on positive definite 2x2 Hermitian matrices and 

B=(l°) , M=(;.*) and / = ( ° I ). (18) 



o -l j > \ z * € J V _1 / 

It was found that in order to obtain a finite limit e — > 0, the normalization constant has to be chosen C e ~ 1/ loge. 
The appearance of this logarithmically diverging term is not surprising since the inverse fermion determinants are 
regularized as 

def 1 ^ + ^70 + z)det- 1 (-D 7llv + M7o + z*) = lim det" 1 ( + ^ D " 7 " + / 7 ° + Z ) . (19) 

Zero eigenvalues of D^j^ do not lead to zero eigenvalues of the operator in the r.h.s. of this equation. Therefore, 
generically its eigenvalues are nonzero real and occur in pairs e ± , and we expect that the e dependence enters 
in a similar way to the mass dependence of the QCD partition function with one bosonic quark at zero topological 
charge. In the microscopic limit this is given by K^ijri) which diverges logarithmically in m. 

To derive the partition functions |Q and (fT7|l we have only used the global symmetries of the underlying QCD 
partition functions. In particular, this implies that any microscopic theory with the same global symmetries and 
spontaneous breaking thereof will have the same zero momentum effective theory. The simplest theory in this class is 
chiral Random Matrix Theory at nonzero chemical potential which is obtained from the QCD partition function by 
replacing the matrix elements Dirac operator by (Gaussian) random numbers. The partition function 117|l was also 
derived explicitly starting from a chiral random matrix model instead of only using symmetry arguments 26] . Another 
advantage of using a random matrix model is that one can easily perform numerical simulations. For example, the 
quenched spectral density was calculated numerically 26] and was found to be in agreement with (|16l) . 

Below we will show that the replica limit of the Toda lattice equation H16|) can also be used to obtain the unquenched 
spectral density using a similar regularization of the partition function with bosonic quarks. In this case we need 

z N f ,n=- 

with Nf 0. This involves an integral over the supergroup Gl(Nf\2). This integral is not known explicitly 
for Nf ^ 0. For that reason we will derive Z„ f ' n ~ 1 directly from the corresponding random matrix model instead 
of the low-energy effective partition function based on Gl(Nf\2). 



III. RANDOM MATRIX MODEL FOR QCD AT NONZERO CHEMICAL POTENTIAL 

In this section we define the random matrix model which we will use to calculate the generating functions. Random 
matrix models for QCD originally Q focused on the quark mass dependence at zero chemical potential. Later 0, 
a random matrix model including the chemical potential successfully explained why quenched lattice QCD at zero 
temperature has a phase transition at a chemical potential of half the pion mass. However, a disadvantage of this 
model is that no eigenvalue representation is known which is required for the use of orthogonal polynomial methods. 
The random matrix model that will be used in the present paper was introduced in [34| • It differs from the model in 



2 Notice that both in the l.h.s. and the r.h.s. the overall factors z and z* are cancelled so that the l.h.s. of (15) gives the spectral density 
without the zero eigenvalues. 
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by a different form of the chemical potential term. Its partition function can be reduced to an eigenvalue representation 
and allows for an explicit solution in terms of orthogonal polynomials |34| . Because this model captures the correct 
global symmetries of the QCD Dirac operator, its low energy limit is also given by 1)14(1 . A related model defined in 
terms of a joint eigenvalue distribution was introduced in |37| . Although this model is not in the universality class of 
QCD partition functions, the complex orthogonal polynomial methods that were developed for the derivation of the 
spectral density are applicable to the random matrix model introduced in [34j . In the section we will comment on 
the relations between the different models. 

The random matrix partition function with Nf quark flavors of mass to/ and n replica pairs of regular and conjugate 
quarks with masses y and z* each is defined by 

f if/ 
2N f ' n ({m f },y,z*;fi) = / d$d* u> g ($)u;g(*) J] det(Z>( M ) + to/ ) 

/=i 

xdet n (2%) + t/)det n (X> t (/z) + ,**), (20) 
where the non-Hermitian Dirac operator is given by 

i$ + /x* 
1<P' + /it 1 



Negative numbers of flavors Nf,n < will denote insertions of bosons given by inverse powers of determinants. To 
allow for the presence of both bosonic and fermionic flavors we will distinguish between Nt> and Nf in later formulas. 
Here $ and \& are complex (N + v) x N matrices with the same Gaussian weight function 

w G (X) = exp( — NTtX^X ) . (22) 

The matrix T>(/i) replaces the QCD Dirac operator plus chemical potential in (01. It has exactly v zero modes 
which identifies v as the absolute value of the topological charge. The model (|2(JII has the same flavor symmetries 
as QCD which are broken by the quark masses and the chemical potential in exactly the same way as in QCD. The 
only difference with the model Q is that the chemical potential is multiplied by a complex matrix \1/ with Gaussian 
distributed matrix elements. In the microscopic limit, where /J?N, to/TV and zN are fixed as N — > oo this partition 
function will match the partition function (|14|) after fixing the scale of the parameters appropriately. The latter is 
uniquely determined in terms of the flavor symmetries and their breaking, and thus universal. As was already done 
in the previous section, the microscopic limit of the partition function will be denoted by Z instead of the notation Z 
which we use both for the QCD partition function and the random matrix theory partition function. Generally, the 
overall normalization of Z and the microscopic limit of Z is different. 

The form of the Gaussian weight wg{&) respects the flavor symmetries, but these symmetries do not exclude traces 
of high er p owers of in the exponent as well as other non-invariant terms. At zero chemical potential it was 
shown [5(1 mi 02, liOi DHL HE iMA H3 * nat partition functions and eigenvalue correlations in the microscopic limit are 
independent of the form of this weight, which, in random matrix theory, is known as universality. The only condition 
for the probability density is that the spectral density near the origin is nonvanishing in the thermodynamic limit. 
In the derivation of the joint eigenvalue density of l|2U|) it is essential that the two weight functions are Gaussian. 
However, we believe that this is only a technical requirement and higher order invariant terms will not alter the 
microscopic limit of the joint probability distribution (for more discussion see section E} - 

In |34| it was shown that the partition function (|20() has an eigenvalue representation after choosing an appropriate 
representation for the matrices <!> and 4" yielding 

. N 

Z^' n {{m f },y,z*^)^{yz*Y n \{m) / Y[d 2 z k V^' n ({ Zi },{z*};a), (23) 

f=i Jc k=i 

where the integration extends over the full complex plane and the joint probability distribution of the eigenvalues is 
given by 

1 N 

^" n ({*},K};a) = — \A N ({zf})\ 2 [J zfc a). (24) 

M fc=i 

The Vandermonde determinant is defined as 

N 

A N ({zf}) = J] (z? - z)), (25) 

i>j=l 
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and the weight function reads 



Nf 



xexp(-ffi-^l^ + z fc * 2 )), (26) 

where if,, is a modified Bessel function. The parameters of the weight function, {mA,y, z* and /i are collectively 
denoted by a. In the quenched case, Nf = n = 0, the weight function will be denoted by w(zk, zjj!; n). We have included 
a normalization factor 1/ fi 2N in eq. 12411 . This way it reduces to the joint eigenvalue density of the chGUE times a 
product of delta-functions in the limit fj, — > 0. We expect that the terms that are nonvanishing in the microscopic 
limit are universal and are required for the derivation of the microscopic spectral density. 

In the asymptotic limit, N\zk\ 2 / n 2 ^> 1, the modified Bessel function in the weight is well approximated by its 
leading order asymptotic expansion. The corresponding joint eigenvalue distribution was introduced in 37] as a model 
for QCD at nonzero chemical potential. Let us look closer at the condition N\zk\ 2 / fJ 2 3> 1- In terms of the half-width, 
i mra ~ /i 2 , of the cloud of eigenvalues and the spacing, A ~ 1/N, of the imaginary part of the eigenvalues it can be 
rewritten as 

|z fe | 2 /(a; max A) > 1. (27) 

For a; max > A this implies that \zk\ ^> A where we expect to find the weight function of the Gaussian Unitary 
Ensemble perturbed by an anti-Hermitian matrix which indeed has a Gaussian weight function (see ,36]). In other 
words, the appearance of the modified Bessel function K v is a direct consequence of the chiral symmetry of the 
problem. 



IV. GENERATING FUNCTIONS FROM RANDOM MATRIX MODELS 

The aim of this section is to derive explicit expressions for the generating functions introduced in section ITU using 
the random matrix model introduced in the previous section. In particular, we will present new expressions for the 
corresponding partition functions with any given number of fermions and one pair of conjugate bosonic quarks, as they 
are needed for the replica limit of the Toda lattice equation. We will find that the partition functions with bosonic 
replicas have to be regularized. This divergence will lead to a factorization into a partition function containing only 
fermions and a purely bosonic one. After taking the microscopic limit, they can be expressed entirely in terms of 
known partition functions 114|) and l|17l) . In the following section we consider partition functions with Nf flavors and 
n fermionic replicas to show that the random matrix model partition functions satisfy the Toda lattice equation 1)150. 



A. Orthogonal polynomials and Cauchy transforms 

As is the case for Hermitian random matrix models, the partition functions 112311 can be expressed in terms of 
orthogonal polynomials, and the integrals to obtain the spectral correlation function can be performed by means of 
the orthogonality relations. In this section we present the polynomials corresponding to the weight function l|26(l for 
Nf = n = (which is denoted by w(z, z*; /x)). 

The orthogonal polynomials are defined as solutions to the following orthogonality relation 

J d 2 z w(z,z*;fi) p k (z) pi(z)* = r k 8 k i , (28) 
c 

where the integral is over the full complex plane with measure d 2 z = dRez dlmz, and r k denotes the squared norms 
of the polynomials. Since the quenched weight w(z, z*; /i) appearing in this relation is real positive the polynomials 
form a complete set with positive squared norms r k ■ The polynomials p k (z) depend only on the variable z and not 
its complex conjugate, as is indicated by the notation. They are also functions of the chemical potential fi, the size 
N of the random matrix, as well as of the topological charge v. We suppress the dependence on /i, N and v. 
In monic normalization, Pk{z) = z 2k + . . ., the solution of l|28|l is given by [34| 
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Since the Vandermonde determinant only contains squared variables only polynomials in those squared variables will 
appear. The same happens for the polynomials on the real line when mapping the Laguerre ensemble from the positive 
real numbers to the full real line. For the weight (|26|l the norms are given by 

TT^ 2 (l+Zi 2 ) 2 ^ fc! (k + V )\ 
T k = N 2k+ V +2 ' ( 30 ) 

Because the Laguerre polynomials L^\z) have real coefficients it holds that Pk(z)* — pu{z*). 
The Cauchy transform of the polynomial Pk(z) is defined by 

h k (y)= [ d 2 z w(z,z*;p.) p k (z)* 1 (31) 

Jc z - y 

We note that h (y) has a nontrivial y-dependence, unlike the polynomial po(y). In contrast to random matrix models 
with eigenvalues on a curve in the complex plane, the integral over the pole in the complex plane is always well 
defined. However, if we take the limit /j, — > of an anti-Hermitian Dirac operator we have to make sure that the 
argument y lies outside the support of eigenvalues. 

The leading order asymptotic large-y behavior of hu[y) is obtained by expanding l/(y 2 — z 2 ) in a geometric series. 
Because of orthogonality all coefficients up to order l/y 2k vanish resulting in the asymptotic behavior 

My) = -^ + o(^)- (32) 

This result is useful for checking the identities derived below. 



B. Partition functions with one flavor 



The polynomials and their Cauchy transform enjoy the following relations to one flavor partition functions. The 
former is related to a single fermionic flavor partition function through the so-called Heine formula, 

-Z'a/ ( z i m) _ / i r 2 » . \ ,. , . 

/-'■<- Zj ) ) = z p N {z) ■ (33) 




„JV f =0 - \ 1JA~ J - 

N 



Here we have written the fermion determinant as an expectation value with respect to the quenched partition function 
at Nf = n = 0. (Vanishing replica (or flavor) indices will be commonly suppressed throughout the following.) For 
/i < 1 the partition function is manifestly real and positive for real masses z, as the orthogonal polynomials are 
Laguerre polynomials of argument —Nz 2 /(1 — /i 2 ) with zeros on the imaginary axis. For fi > 1 the zeros are located 
on the real axis. Remarkably, the Yang-Lee zeros [5^| of the model in behave quite different from the model l|20() 
discussed in this paper. The polynomials can also be interpreted as characteristic polynomials. 

The relation (|33[1 can also be written down for Pi^n{z), expressing it as an average over I variables. However, in 
order to be precise we would have to chose the weight function in (|26|l to be TV-independent as in [3^, |3^, ^3 • Since 
we are mainly interested in taking the large- N limit we prefer to explicitly keep the N or volume dependence inside 
the weight. For I N we can safely ignore this technical subtlety. 

Turning to the Cauchy transform of pn(z), it is given by the partition function with one bosonic quark flavor 

Z " \ ^ / Nf =n=0 ^ 

The same remarks about positivity and the A^-dependence of the weight function made before also apply here. 

Our program is to express spectral correlation functions in terms of partition functions. This has the advantage 
that the microscopic limit is given by a group integral based on the symmetries of the partition functions. To achieve 
our goal we also need all possible two flavor partition functions. At zero chemical potential the pro gram of expressing 
spectral correlation functions in terms of multiflavor partition functions was carried out in 



:nc prog 

mil. 
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C. Partition function with two flavors 



It is known from Hermitian random matrix models |6l| that the expectation value of two characteristic polynomials 
is proportional to the kernel of the orthogonal polynomials i|28|) . As was shown in j3£|, for an arbitrary weight 
function, similar results can be derived for non-Hermitian random matrix models. Since we are now dealing with 
non-Hermitian matrices we have to distinguish between a characteristic polynomial (or mass term) and its Hermitian 
conjugate. For a single pair of one fcrmionic quark with mass z and one conjugate fermionic quark (with complex 
conjugated eigenvalues) with mass u* we have the relation 38] 



Z " iyj^ = ( (zurf[(z 2 -zf)(u 2 -z* 2 ) ) = r N (zuYK N+1 (z,u*), (35) 

Z N \ J=l / Nf =n=Q 

N 1 

K N+1 (z,u*) = 22—Pj(z)pj(u*). (36) 



3=0 J 

Here we have defined the bare kernel K^(z,u*) which differs from the full kernel by the square root of the weight 
function at each argument. The full kernel appears in the computation of correlation functions from orthogonal 
polynomials (see section Ivj) . In particular, for z — u we obtain the result for one pair of fermionic replicas n = 1 

gr^iM) = / |z|2 ,-Q |z2 _ Zj2|2 \ = rN \ 7rKN+i[ZiZl (37) 

Z N \ 3=1 / Nf =n=0 

This quantity is manifestly real and positive. On the other hand, if we compute the partition function with two flavors 
with the same eigenvalues, we obtain |38| 



2 ( X1 ,X 2 ;H) _ / I -TT, 2 2w 2 2\\ {xiX^Y 



N,=n=0 



Pn(xi) pn+i(xi) 
Pn(x 2 ) Pn+i(x 2 ) 



(38) 



For Hermitian random matrix models the two results (|35|) and l|38ll are related through the Christoffel-Darboux 
formula, which in general does not hold in the complex plane (see also | Appendix A"| ). 

We now turn to the most important two-flavor partition function, containing a pair of a bosonic quark and its 
conjugate. In the replica limit of the Toda lattice equation this partition function appears when computing the 
quenched density (see i|16|l ). In the microscopic limit this partition function was calculated in |2fij |. both by starting 
from a random matrix model, and by integrating explicitly over the Goldstone manifold in the chiral Lagrangian i|17[l . 
It turned out that this partition function is singular and has to be regularized according to the prescription given 
in (|19|l . The partition function is obtained in the limit of vanishing regulator, e — > 0, where it diverges as ~ log(e). 
This diverging constant can be absorbed in the normalization of the partition function. Below we will establish the 
nature of this singularity in the orthogonal polynomial approach. In fact, it does not enter in the calculation of the 
spectral density for the random matrix model l|20|l by means of complex orthogonal polynomials |34|. However, we 
will show that it occurs even at finite N in the random matrix model partition function with one bosonic quark and 
one conjugate bosonic quark of the same mass. 

Generalizing the results in |39J to the chiral case we find that the following relation holds 

zr° = \ {yx) -*? 9 ) / = An ~ 2{x (39) 

N-2 ^ 

A N ^ 2 (x*,y) = -Q{x*,y;n) + J2— h j( x *) h j(y)> ( 40 ) 



r 

3=0 • 



where 



( ~ )lr = / d 2 zw(z lZ *-^ ^ 2 _ x ^ z2 _^ . (41) 



c 



Equation (|40|l defines a bare kernel containing only Cauchy transforms. For different masses x ^ y both the sum over 
Cauchy transforms and the function Q(x* fi) are finite. At equal arguments, x — y, the sum remains finite but 
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the integral in l|41|l is logarithmically divergent, just like the divergence encountered in i|17|) . In order to regularize 
this divergence we cut out a circle C(x, e) of radius e around the pole at z — x and one around the pole at z — —x. 
To evaluate the singular part of the regularized Q e (x* , x; fi) we expand the weight function around the singular point 
and keep only the leading part. For \x\ 3> e we find for the pole at z = x 



Q e (x,x*; n) 



w(x, X* 
w(x, X* 
w(x, X* 



dzdz* 



C-C(x 



C-C(x 



-2% {z* 2 
dzdz* 1 
-2i 2x* 



x* 2 )(z 2 — x 2 ) 

log((z* -x*)/(z* +x*)) 



dz 1 log((z* -x*)/(z* +x*)) 



ttw(x, x* 



2xx* 



2i 2x* 
loge + O(e ). 



dC(x,e) 

M) 



(42) 



We find the same contribution from the pole at z = —x. Thus we have obtained the following interesting relation for 
the singular part of the bosonic partition function 



7TL — — 1 



7 N f =0 



(M) 



irw(x, x*; jj,) 
~2\x\ 2 + 2v rjv-i 



loge + 0(e°) 



(43) 



We are not aware that such a relation between the weight function (|2(j|l of complex orthogonal polynomials and 
the bosonic partition function at a degenerate mass pair was previously known. It translates easily to non-chiral 
non-Hermitian random matrix models. 

Relation 1)43(1 shows that the singularity that was found in [2|| for the microscopic limit of Z^ J ~~ 1 (x 1 x*; /i) even 
persists for finite size random matrices and can be expressed simply in terms of the weight function. Only weak 
assumptions have been made in the derivation of (|43|l . and its validity extends to a wide class of weight functions. For 
x r~j 0(e) the above derivation does not apply. Instead we find that lim^^o Qe{ x , x*; /i) ~ 1/e 2 - As we will explain 
in the next section, the logarithmic singularity we have found for the two-point function l)43|l will also occur in more 
general partition functions with at least one pair of conjugate bosonic quarks. 

The last two-point function we need below is the partition function with one fermion and one boson. This partition 
function reads 



N 



l-.N h 



7 N f =0 
J N 




-j-j- (x 2 z 2 ) 




( y 2 - x 2 ) - M N -x[x,y) 



— Pj(x)hj{y) . 



(44) 



(45) 



Here, we have defined the bare kernel A/jv-i(^) y) consisting of orthogonal polynomials and Cauchy transforms. This 
partition function is nonsingular for any values of the arguments and correctly normalizes to unity at equal arguments 
y = x. 

The remaining two flavor partition functions combining (conjugate) fermionic quarks and (conjugate) bosonic quarks 
are given in |Appendix A| 



D. Partition Functions with an arbitrary number of flavors and one conjugate flavor 

General expressions for the expectation value of ratios of characteristic polynomials valid for an arbitrary complex 
weight function were given in (3^, H^, . We will apply these results to partition functions with Nf flavors plus one 
pair of regular and conjugate fermionic (n = +1) or bosonic (n = —1) flavors which both enter in the expression for 
the spectral density l|16fl after taking the replica limit n — > 0. 

The partition function with Nf fermionic quarks with masses m/ is given by |38| 



N t(S \ \ I Nf N \ U N} m v 

*N \ f=l 3 = 1 I N,=n=0 /U fS ' ~ ~ 



11 



Next we add a single pair of a fermion and its conjugate. This partition function is given by [38 



Z 



Nf,n=+1 
N 



({m f },z,z*;n) _ 



T N f =0 
J N 



N f N 

nn-/)n 

/=1 3=1 



N, 



,]| 2 II(-/-^ 
/=1 




N f =n=0 



{-l) N 'r N \z\ 2 "Ufii m f 
A Nf+1 ({m 2 f },z 2 ) 



PN+i(mi) 



PN+i{m Ns ) Pn+i(z) 



PN+N f (mi) . . . p N+Nf (m Nf ) PN+N f (z) 
K N+ i(mi,z*) ... K N+1 (m Nf ,z*) K N+1 (z,z*) 



(47) 



where Km+i{x : y) is the bare kernel introduced in QHfiJI. In the replica limit of the Toda lattice equation, the partition 
function with one pair of conjugate bosons instead of conjugate fermions enters as well. It is given by |8fll | (see also 
|Appcndix A| | 



7 N f ,n=-l 
~*N 



({m f }\y,x*;n) 



n 



Nf v N 
f=l m f 



7 N f =0 
"A' 



n 



Ufii m f( m 'f - v' 2 ) 

"rjv_i A Nf ({m 2 f })(yx*) v 



\ {yx*Y fJi \ (z]~y 2 )(z 
PN-i(mi) ...p N -i(mN s ) 




(48) 



N f =0 



PN+N f -2{mi) 

J\fN+N f -2(mi,y) 



PN+Nf-2(mN f ) h,N+N f -2{x*) 

J^N+N 1 -2{m Nl ,y) AN+N f -i{x*,y) 



If we take the limit of degenerate bosonic masses, y — > x, the function Q(x* ,y; /i) inside the matrix element 
AN+N f ~2{%* , y) becomes singular, while all other matrix elements remain finite. Expanding the determinant with 
respect to the last row this singular part ~ Q e (x*,x; fi) gives the dominant contribution, 



7 N S =Q 
J N 



det 



r N _i A Nf ({m 2 f })\x\ 2,/ i<k,l<N f 



PN+k-2{mi) } A N +N f -2{x* ,x) . (49) 



Using the relations (|40|l and (|46|l we thus arrive at the following factorized result, 



lim Z ^ - 1 ({ m f}\ x > x *>^ 



7 N f =Q 



/=i 



x 2 ) 



■ 7 N f =Q 
-°N-l 



lim 



z^—Hx.x*-^) 

„N f =0 
N 



(50) 



Both sides are regularized by replacing Q(x,x*;n) by Q e (x, x*; fi) as discussed above |@2J. This is the main result 
of this section. Let us stress that the factorization we obtained is not a consequence of the suppression of non- 
planar diagrams in the large N limit. In (|50|l it occurs at finite- AT and is strictly due to the singularity for coinciding 
arguments. This completes the derivation of all generating functions required for the derivation of the spectral density 
from the replica limit of the Toda lattice equation (|16fl . 

If we were to compute not only the spectral density but also higher order correlation functions from the Toda lattice 
approach (see e.g. in j2||), we would have to add more pairs of replicated bosons, each with different masses yj and 

x* . Taking the degenerate limit yj — * Xj would lead to the same type of factorization into ~ FT ■ Z^ ~ 1 (xj ,x*;fi) 
and a remaining fermionic partition function. 



V. TODA LATTICE EQUATION FROM ORTHOGONAL POLYNOMIALS 

The purpose of this section is threefold. First we show that the generating functions with n pairs of fermionic 
replicas satisfy the Toda lattice equation. Since the microscopic limit of the random matrix partition function, Z^ f ' n , 
is given by the effective partition function (|14|) 2(| |3^| , this result also shows that the Toda lattice equation (|15fl for 
Nt = |2g can be extended to arbitrary Nf. Second, we show that the spectral density obtained from the replica 
limit of the Toda Lattice equation agrees with the result computed from orthogonal polynomials |34| . This extends 
the agreement between the two approaches to cases where the weight in the partition function is no longer positive 
definite. We close this section with some remarks about universality. 
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A. The spectral density from the Toda lattice equation 



We start by showing that the Toda lattice equation l|15fl can be modified to hold for random matrix model partition 
functions at finite N. The starting point of the evaluation is the expression for the partition function with Nf + n 
ordinary flavors and n conjugate fermionic flavors generalizing (|47|l . While this result is given in |3S| for non-degenerate 
masses we have to take the limit where the replicated masses become degenerate. This leads to successive derivatives 
in these variables of the polynomials and of the kernels l|3()[) . For simplicity we only display the result for Nf = 1, 



N 



1 (to, z, z*; fi) fn v 



\1nv T-iN+n-1 _ 
1 lj—N 



r N f =0 



(m 2 - z 2)n| z |n(n-l) 



(51) 



p N+n (m) PN+n(z) S z p N+n (z) 

K N+n (m,z*) K N+n (z,z*) 5 z K N+n (z,z*) 
8 z *K N+n (m,z*) 8 z *K N+n (z,z*) S z *S z K N+n (z, z*) 

8 n z ^ l K N+n (m,z*) 8 n ^ l K N+n (z,z*) 5 n ^H z K N+n {z,z*) 
omitting constant factors. Here the derivatives are defined as 



S^Kk+^z*) 



8^ l 8^ l K N+n {z,z*) 



8 Z = z—, 
az 



._d_ 

dz* 



(52) 



We observe that this expression depends on the matrix size explicitly through the combination N + n and implicitly 
through factors that originate from the iV-dependence in the exponent of the weight function. Using the Sylvester 
identity as in 26] while keeping only the N dependence that enters asJV+nwe easily derive the Toda lattice equation, 



5 Z 8 X * log Z% f 1,n 



(to, z, z*; fi) ~ n(zz*)' 



■ 2 n-i {m,z,z 



Im)Z 



N+l 



-1/ 

[m, 



://) 



N f = l,n 
N 



(m,z,z*;fi)f 



(53) 



again ignoring constant factors. This equation is valid at finite N provided that the factor N that appears explicitly 
in the exponent of the weight function is not changed. The subscript of the partition function thus denotes the total 
number of complex integration variables. For large N this distinction can be ignored. It is not difficult to generalize 
the Toda lattice equation to arbitrary Nf. The result is 



8 Z 8 Z * logZx S ' n {{mf},z,z*;fi) ~ n{zz*) 



2 2 N 



N f ,n+1 
1 



(W}> 



tj)Z^ 1 ({m f }, 



(54) 



We still have to fix the overall normalization constant of the partition functions which in principle depends on N, Nf , 
as well as on the replica index n. Since the l.h.s. of (|54|) is linear in n we have to choose the normalization constants 
in the r.h.s. to retain this rt-dependence. With the appropriately adjusted constants the random matrix model 
generating functions satisfy the equality JTSJl. 

Since the microscopic limit of the random matrix partition functions Z^ f ' n ({mf}, z, z*; /i) is given by the group 
integrals (|14l) . we conclude that l|15l) is satisfied for any number of flavors. Explicit expressions for these partition 
functions in terms of Bessel functions will be given in section IVT1 where we also give results for the spectral densities. 

Before doing so let us step back to discuss what we have achieved for the spectral density using the replica limit 
of the Toda lattice equation Combining equations (|15|) and (|T2*|l with the factorization (|5H|l and canceling one 

of the partition functions in the denominator of the Toda lattice equation, we arrive at the following result for the 
spectral density, 



P» f (z,z*,{mf};n) = - 



7r log e 2 



Ifc' 



7 N f ,n=l 
■'N-l 



({mf},z,z*;fx) 



2 N f ({mf}; fi) 



(55) 



where we have explicitly displayed the divergent normalization factor —{tt/2) loge which cancels the divergent factor 
that was obtained in l|43() . The v dependence of the spectral density is contained implicitly in the partition functions. 
This is the central result of the Toda lattice approach applied to the random matrix model (|20|l . The spectral density 
factorizes into QCD like partition functions which, at low energy, can be expressed as group integrals of the form 1141) 
and ifT7jl . 
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Let us discuss some properties of the spectral density. First, it is symmetric under z — > — z (together with 
z* — > —z*), because the partition function Z^ , ' n ~ 1 contains only squared variables and Z^[~ _1 is even as well. This 
is a consequence of the axial symmetry of the partition function which requires that the nonzero eigenvalues occur in 
pairs of opposite sign. Inserting the result Ij47(l for Z N f ' n ~ l we see that the prefactor Il/=ii( m / — z2 ) IS canceled from 
the Vandermonde (for an explicit expression in terms of determinants sec (15711 below). Replacing the polynomials 
and kernels in Z N f ' n ~ by their expressions in terms of partition functions we notice one further symmetry (for real 
quark masses), 

p v f (z,z*,{mf};fi)* = p v f (z*,z,{mf};fi), (56) 

that is taking the complex conjugate and at the same time exchanging z and z* . These symmetries follow directly if we 
write the spectral density as an integral over the joint eigenvalue distribution. Thus the real part is symmetric about 
the real axis in the complex eigenvalue plane while the imaginary part is anti-symmetric. Therefore, the eigenvalue 
density on the real axis is real. However, away from the real axis it is, in general, neither real nor positive and 
therefore does not define a probability density. The lack of reality directly relates to the asymmetric way the variables 
z and z* enter in the polynomials and kernel respectively. The conjugation symmetry is related to the invariance of 
the partition function under changing the sign of /x, and is valid for QCD beyond the ergodic domain J2J as well. 



B. Equivalence of the Toda lattice and orthogonal polynomials approach 



In order to compare the result (|55|l from the replica limit of the Toda lattice equation to the result from orthogonal 
polynomials [34| it is instructive to rewrite l|55(l in terms of the weight function, polynomials and kernels. Reinserting 
(g3J), (g3 and lO wc obtain 



Pv' ' {z,z*;{rrif},p) = w(z,z*;p) 



pjv(mi) 



p N {m Nf ) 



Pn(z) 



PN+N f -i(mi) . . . PN+N f -i(m Nf ) p N+Nf -i(z) 
K N (mi,z*) ... K N (m Nf ,z*) K N {z,z*) 



deti<k,i<N f [pN+k-i{mi) } 



(57) 



again after normalizing appropriately. As we will explain next, this is precisely the result that was obtained in [34| . 

It is well-known from random matrix theory that the spectral density can be written in terms of a kernel of 
orthogonal polynomials [62(. This also holds for complex weight functions. However, if the weight function is not a 
symmetric function of z and z* , we have to use bi-orthogonal polynomials in the complex plane defined by [39ll63|. 



Sz v£<%Xi{m,},ri p[ N! \z) ^V) = r Wff„, 



(58) 



assuming the pseudo norms to be nonvanishing. (Bi-orthogonal polynomials in general no longer form a scalar 

product with positive definite norms.) In contrast to orthogonal polynomials on the real line, the Pk{z) and qk{z) are 
in general different polynomials. Defining their kernel by 



i N ~ 1 1 

K. { x f) (z,u*) = [w^ f ' (z,z*;{m f },fi)w^ f '°(u,u*;{m f },fi)}^ ^ ~JWJ) p j ( u *^ 

j=o r 3 

the following expression for the correlator of k complex eigenvalues holds [6^ | 



pv f {{zi,zt}i=i 1 ... lk ;{m f };ti) 
In particular, for the spectral density we obtain 



det 

l<j,l<k 



N-l 



Pv f (z,z*;{m f };n) = Wv f '°(z,z*;{m f };n) ^ — ^-y p^ Nf) (z)q^ N,) (z*) 



3=0 >j 



(59) 



(60) 



(61) 



In [34| it was shown that the kernel for Nf flavors can be expressed as a determinant of kernels for zero flavors which 
immediately leads to the result l|57|) . 
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This establishes the equivalence between the Toda lattice approach and the bi-orthogonal polynomial approach at 
the level of the spectral density. Starting from the observation that also in the case of generating functions for multi- 
point correlation functions the partition function Z^~^ factors out, we are confident that the equivalence 
between the Toda lattice approach and the bi-orthogonal polynomial approach can be established in that case as well. 

C. Universality 

After having computed the spectral density from the replica limit of the Toda lattice equation and compared to that 
from the orthogonal polynomial approach, we would like to use the insight gained to discuss the issue of universality. 
First of all, we can distinguish between two different ways to approach universality. The first way is based on the 
observation that different theories with the same pattern of spontaneous symmetry breaking and a mass gap have the 
same low-energy limit. For example, gauge field theories and random matrix theories with the same global symmetries 
are described by the same effective partition function in the microscopic scaling limit. For QCD with three colors 
and fundamental quarks at nonzero chemical potential this universal partition function is given by the group integral 
1)14(1. In QCD, a mass gap exists because of confinement, and in Random Matrix Theory, a mass gap appears in the 
large N limit. Therefore, if we change the action without affecting the global symmetries and the existence of a mass 
gap, the microscopic scaling limit will remain the same. A second way to understand universality is based on a direct 
calculation of correlation functions for a deformed probability distribution which does not necessarily have the unitary 
invariance of the Gaussian model. In the case of deformations that respect the unitary invariance, universality then 
follows from the asymptotic behavior of the modified ortho gona l polynomials. In particular, for chiral random matrix 
theory at /i = this approach has been very successful [5lL l55| . 

Let us also comment on the role of the weight function in non-Hermitian random matrix theories. From the replica 
limit of the Toda lattice equation it is clear that the eigenvalue density will be proportional to the partition function 
with one pair of conjugate bosonic flavors which, according to this work, is proportional to the weight function. In 
general, this weight function contains both universal and non-universal parts, and in the microscopic scaling limit 
only the universal parts survive. If we allow for deformations of the Gaussian weight by higher order powers the 
resulting weight function, and therefore the spectral density, will be modified. It is only in the microscopic scaling 
limit that the non-universal parts of the weight function disappear and a universal spectral density will be recovered. 
In Hcrmitian matrix models the weight function at finite N is also non-universal but reduces to a universal result 
in the microscopic scaling limit. Non-Hermitian random matrix models for QCD at fi ^ represent a continuous 
one parameter deformation of the Hcrmitian model at [i = 0, governed by ^jl 2 N . In this case, the weight function 
contains an additional universal piece in the microscopic scaling limit that develops into a delta function ~ (5(Im(z)) 
for [i 2 N — ► 0. The additional universal piece occurs because an eigenvalue representation can only be obtained 
after a nontrivial integration over the similarity transformations that diagonalize the non-Hermitian matrices. In the 
Hermitian case it is trivial to obtain the eigenvalue representation, and no additional universal piece is generated. For 
technical reasons we have not been able to integrate out the similarity transformations when higher powers are added 
to the Gaussian probability distribution of the model (|20|) . We mention in passing that for non-chiral non-Hermitian 
random matrix models a universality proof by deforming the random matrix potential could be given |64| . 

VI. RESULTS AND EXPLICIT EXAMPLES 
A. The partition functions in the microscopic limit 

Having obtained analytical expressions for the partition functions of the random matrix model l|2(JII at finite N we 
now consider the large N limit. This limit will be taken according to the counting scheme where m/N and /i 2 N are 
kept fixed as N — > oo and is referred to as the microscopic limit. The rescaling of the non-Hcrmiticity parameter 
/i 2 with N was first introduced in |65j| as the concept of weak non-Hermiticity. As discussed in the previous section, 
in the microscopic limit the random matrix partition function is given by an integral over the Goldstone manifold, 
i.e. it is uniquely determined by the spontaneous breaking of the flavor symmetries. The random matrix model can 
be viewed as an alternative way to calculate the integral over the Goldstone manifold, which in the case of nonzero 
chemical potential and both bosonic and fermionic quarks has not been accomplished in any other way. 

The dependence of the quark mass and the chemical potential in the integral over the Goldstone manifold is only 
through the combinations mlW and /i 2 F 2 V, cf. (|14fl and (|17|l |14j. The dimensionful scales can be recovered from 
the random matrix model by making the replacements 



2mN -> mVE 



(62) 
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2/j 2 N -» n 2 F%V. 

To make contact with the physical content of the theory we will express our results in terms of the dimensionful scales 
in the remainder of this section. 

In the case we have only either fermionic or bosonic quarks and no conjugate quarks, the Goldstone bosons are 
not charged with respect to the chemical potential, and we should find the microscopic partition functions at zero 
chemical potential. One easily observes that the effective partition function (|14fl does not depend on the chemical 
potential in this case. In the microscopic limit of (|33|> we thus find the partition function (|14f> with one fermionic 
flavor [Hililil 

Z? f=1 (m;n) = J v (mSV). (63) 

An overall scaling factor e - '' F -* v has been removed. This unphysical factor can be removed in the random matrix 
model by scaling the Gaussian potential in the weight (1221) by a factor of 1 — /i 2 34] . We did not include an explicit 
scale factor in the weight function, but will assume this procedure has been done in order to arrive at the correct 
expressions in the microscopic limit. 

When there is one pair of conjugate fermionic quarks, the microscopic partition function is given by the result 
obtained by performing the integral l|14|) over the Goldstone manifold |26j ] 

Z" =1 (y,2*;/z) = e 2 ^ F i v / dtte- 2 ^ F * vt2 I„(yY;Vt)l4z*ZVt). (64) 

Jo 

For the universal partition function with a pair of conjugate bosonic quarks we find using (14311 

where we have divided out the divergent factor — 7rlog(e)/2. This is again in agreement with the result [2^| obtained 
by an explicit integration over the Goldstone manifold l|17jl. 

All other microscopic partition functions required to calculate the spectral density can be obtained from the relation 

N f 

Z^' n —\{m f }, z, z*; M ) = J] K ~ s 2 )^' ({"»/}; **;/*)• (66) 



B. The microscopic spectral density 



In this section we write out several different examples for the eigenvalue density of the QCD Dirac operator in 
the microscopic limit. We emphasize that these results are nonperturbative analytical predictions for the QCD Dirac 
spectrum at nonzero chemical potential which should be reproduced by lattice QCD simulations. In order to get 
more compact expressions we will express our results in terms of the partition function (|64|l . The normalization of 
the quenched density is chosen such that there are T,V/n eigenvalues per unit length along the imaginary axis. The 
density of the projection of the eigenvalues on the imaginary axis is therefore equal to the eigenvalue density at /i = 0. 
The other densities reduce to the quenched density when the quark masses are taken to infinity (rrif'EV ^S> fi 2 F 2 V) 
and we use this limit to fix the normalization. 

The general expression for the spectral density is given by 

* N f 7 N f> n=1 (( \ * \ 

p?f(z,z*,{m f } W ) = fi- JJK - *)Zr-Hz,fw) " d,?',* ^ . (67) 

z f =i z » 

where the microscopic limit of Z n N — 1 {z 1 z*-^) is given in the previous subsection. The microscopic limit of 

„N fn ! \ det[(x k d Xh ) l I v (x k )]k,i=o,...,N f -i , CQ . 
Z v t {{m f y,n) = -^-j- '- — , (68) 
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Re[z}ZV 25 

FIG. 1: The quenched eigenvalue density with a chemical potential fiF^y/V = 0.1 (upper) and fj,F w VV = 2.5 (lower). The 
density is real and positive and follows in the other quadrants by reflection on the axes. In the upper figure fiF^y/V <C 1, 
and one can still see the individual eigenvalue distributions. In the lower figure the repulsion of eigenvalues from the origin is 
apparent. The width of the strip is on the order Re[z]E ~ 2fi 2 F%. The normalization is such that the width times the hight is 
independent of /i. 

where Xk = wikVY,. Note that with only regular fermionic quarks, \i has no effect on the microscopic partition 
function. The partition function Z^ 1 ' n ~ 1 ({mf}, z, z*; fi) is given in terms of the polynomials pN{m) and the bare 
kernel Kn+i{x, y) in (|47|l . The microscopic limit of this partition function is obtained by expressing the polynomials 
and the kernel in terms of partition functions that are known in the microscopic limit. The polynomials pAr(m) can 
be expressed as the partition functions with one fermionic flavor (see which in the microscopic limit does not 

depend on [i and is given by I l/ {mVYi) in l|63[) . In (|37|l the kernel Kn + i(x 7 y) was shown to be equal to Z™ =1 which 
is given in l|04|) . In the remaining sections we will set v = for simplicity. 
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1. The quenched spectral density 

The quenched microscopic eigenvalue density of the Dirac operator at nonzero baryon chemical potential is given 
byH 

>~-°<- •»> - ^-w'W*, «^ m 

This result has been checked by direct numerical simulations [2(| within a random matrix model, and successfully 
compared to results from lattice QCD |^. Plots of the quenched spectral density are given in Fig. ^ 

By using the leading order asymptotic forms of the Bessel functions, one easily shows that in the limit /j, 2 F 2 V ^S> 1 
the support of the quenched spectrum is inside the strip [2H, |49| 

l Re ^)l £ < l (70) 

with a density that is equal to 

In Fig. [l\we observe that this plateau is already present for [iF^^/V — 2.5. This step in the spectral density is 
associated with a phase transition of the generating function. The theory with a pair of conjugate fermionic flavors of 
mass x, described by Z n= i(x, x; /i), is in the normal phase for xTi/F 2 = rn 2 /2 > 2/ji 2 and in a Bose condensed phase 
for xT,/F 2 = m 2 ./2 < 2/z 2 For the x-integrated quenched spectral density we find (z = x + iy) 

dx p^lo°(z 7 z*;n) = (72) 

consistent with our choice of normalization. 



2. The spectral density for one flavor 

The microscopic eigenvalue density with one dynamical quark of mass m and with baryon chemical potential \i is 
given by [34| 



P N ^(z z* m-u) l£*Yle-*f&K f^^l --^^ 
Pu=0 [z,z ,m,fi) - 2 ^ 2p 2 e K Q y J 



e 



I (mEV) Io{zSV) 
Z™=o (m, z* ; fi) Z^zl (z,z*; (J,) 



J (m£V) 



I [mHV)Z^[z,z*;n) 



This result is shown in Fig. Hfor p,F^W = 0.1 and mVT, = 10, and in Fig. [J] for ^F^W = 2.5 and mVE = 5. It is 
our first example of a spectral density that is complex due to the sign problem. 

Since the unquenched spectral density is proportional to the quenched spectral density its support cannot go beyond 
the support of the quenched spectral density. In particular for fi 2 F 2 V ^> 1 the unquenched spectral density is also 
confined inside the strip given by l|70l) . As we can observe from Fig. 0] the asymptotic behavior is much more 
complicated than in the quenched case. An asymptotic form of the integral l|64H was derived in . For p?F 2 V ^> 1 
and (x + m)T, / (4fj, 2 F 2 ) < 1, it is well approximated by 

1 roc 

dtte- 2 ^ F * Vi? I ((x - iy)EVt)I (mEVt) « / dtte- 2 ^ F * Vi? I ((x - iy)ZVt)I (mEVt) 

Jo 

* ^ Io A ./> . (74) 



Afi 2 F 2 V V V-^l 

This approximation breaks down if (x + m)E / '(4^ 2 F 2 ) > 1 because then the saddle point of the t-integration will 
be outside the interval [0, 1]. The exponential functions in the asymptotic behavior of the Bessel functions no longer 
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FIG. 2: Real and imaginary parts of the eigenvalue density for one flavor of mass mVY, = 10 and chemical potential 
IJ,F n VV =0.1. With these values the real part of the density is quite similar to the quenched density. Note that the scale of 
the imaginary part is considerably smaller than that for the real part. The imaginary part is odd in Re[z] as well as in Im[z]. 

compensate each other in the additional contribution to the unquenched spectra density. Instead of a plateau for 
Nf — we find an oscillatory contribution with an amplitude that increases exponentially with the volume. For 
the total spectral density a plateau is still visible in the region where the quenched contribution dominates. As an 
illustration we show in Fig. 0]the real part of the eigenvalue density for n 2 F%V — 100 and mVY* = 100. The period 
of the oscillations is of the order of the level spacing at [i = whereas the amplitude is of the order exp(10) for our 
choice of parameters. 





FIG. 3: Real and imaginary parts of the eigenvalue density for one flavor of mass mVT, — 5 and chemical potential [iF n y/V = 
2.5. With these values the real part of the density deviates substantially from the quenched density. Note that the real part of 
the density changes sign at z = m. In this case the scale of the imaginary part is comparable to that of the real part. 



3. The spectral density for two flavors 
For two flavors of mass mi and m 2 at nonzero baryon chemical potential the density of eigenvalues is given by 



|z| 2 E 4 F 3 



\z^V_ 

VF2 



Jo(miEV) Io(m 2 ZV) I Q {zYV) 

mi2Vh{mi£V) m 2 EVIi{m 2 T,V) z?VIi{zYV) 
Z"=a(nii,z*;iJ,) Z™=£(m 2 , z*; n) Z^(z,z*; fj,) 



/o(miEV) miEyJi(miEV) 
/o(m 2 Ey) m 2 £Vh(m 2 ZV) 



(75) 



Nf=l 

Re[p^j (z,z*,m;ii)] 

WV 1 



.001 



-0.001 
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FIG. 4: The real part of the eigenvalue density for one flavor of mass mT,V = 100 and chemical potential (iF^VV = 10. The 
severe sign problem manifest itself in the strongly oscillating region of the eigenvalue density which has amplitudes of the order 
exp(10). The full size of the peaks has been clipped for better illustration. 



In Fig. [S]we show the real and imaginary parts of the spectral density for iiF^^/V = 2.5 and m\YN = m^TiV = 5. 
The spectral density is neither real nor positive. 

4- The spectral density with one pair of conjugate flavors 

The microscopic spectral density of the Dirac operator in QCD with one pair of conjugate quarks (with mass m 
and m*) is found to be (see (|57J) and pBSf) 



71 — 1/ * \ 

P v =o{z,m,m ;n) 



|z| 2 E 4 l/ 3 



27TAt 2 F2 

p v U (z;w 



z\ 2 Y?V 



(z z + z* Z )^ Z V 

8 ,,' 2 I?' 2 



%u=o (m,m*; n) Z^ = l (m,z*;n) 



Z r J=o (z,m*; ji)Z™=l (to, z* ; /x) 
Z^(m,m*;fj,)Z^(z,z*;fi) 



(76) 



From this expression it follows that the density is real, positive, and has a zero at z = to. This is illustrated by Fig. 
El where we plot the spectral density for mVT, = 5 and /iF^VV = 2.5. 

Since the quark mass, to, is real, the determinant of the conjugate quark is identical to the determinant of a quark 
with a chemical potential of the same magnitude but with opposite sign, cf. That is, the partition function 

in this case is equivalent to a theory with two degenerate flavors and nonzero isospin chemical potential, /ij. This 
version of QCD does not have a sign problem [6^ as has been explored by lattice QCD simulations |46| . 

5. The density with one ordinary flavor and one pair of conjugate flavors 

As a final example we write out the microscopic eigenvalue density in QCD with 3 light flavors with real masses 
m u = JTid = tti and m s and chemical potentials \ii = fx s = /j. This theory has a sign problem due to the extra quark 
that isn't matched by a conjugate quark. The phase diagram as a function of \ii and /i s was determined in [7f]j| . The 
eigenvalue density is 



N t =l,n=l, * s. 

P v =o K z i z ,m s ,m,m;fi) 



-2fj, 2 F^V 



VF 2 
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FIG. 5: Real and imaginary parts of the spectral density with two flavors of mass miVE = m-zVYj = 5 and chemical potential 
fiFnVV — 2.5. The density is neither real nor positive and the fluctuations are an order of magnitude larger than the quenched 
density at the same value of fi. The spectral density has a zero at z — mi. The imaginary part of the spectral density is 
anti-symmetric about the real and imaginary axis. 



J (mEV) I (m s ZV) h(z 
Z^(m,z*;n) Z^(m s ,z* ]fJ ,) Z^(z 


,171', fi) 




I {mZV) I (m s ZV) 
Z"=o (m,m;n) Z™=1 {m s ,m; /i) 





(77) 



For z = m the two top rows are identical and the first and third column are identical in the 3x3 determinant. Hence 
the density is zero at z — m but does not change sign. While for z — m s the first and second column in the 3x3 
determinant are identical, so the density changes sign at z = m s . A plot of this density is shown in Fig. [7] 




FIG. 6: The spectral density with one pair of conjugate flavors of mass miVT, = U12VE — 5 and chemical potential fiF n y/V = 
2.5. It is real and positive with a zero at z — mi. Since the masses are real and equal this corresponds to a nonzero isospin 
chemical potential. 

N f = l,n=l 

Re t/V=0 (z,m,m,m s ;n)} 
WV 2 




FIG. 7: The real part of the density in a theory with 3 quarks. The masses are m u V"E = m^VY, = 5 and m a VT, = 10. The 
chemical potentials are fj, s = fii = (jl = 2.5/(F 7I \/V). Note that the density bounces off the real axis at z = m u and changes 
sign at z = m a as expected. 

VII. CONCLUSIONS 

We have analyzed the spectrum of the QCD Dirac operator at nonzero baryon chemical potential. In the microscopic 
limit this spectrum is uniquely determined by the global symmetries of the QCD partition function. This is true both 
for the quenched and the unquenched theory. In both cases the spectral density of the Dirac operator can only be 
obtained after the introduction of a complex conjugate pair of valence quarks resulting in a nontrivial baryon charge 
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matrix and Goldstone modes with nonzero baryon number. The microscopic limit of the generating function for the 
QCD Dirac spectrum can therefore be represented in two different ways. First, as an integral over the Goldstone 
manifold which is determined by the pattern of symmetry breaking. Second, as the large- N limit of a partition function 
of an ensemble of non-Hcrmitian matrices with the symmetries of the QCD partition function. In recent work, the 
Dirac spectrum was derived directly from the random matrix model by means of complex orthogonal polynomials. 
In this paper we have obtained the Dirac spectrum from the replica limit of the Toda lattice equation. This explains 
that the spectral density factorizes into the product of a fermionic partition function and the partition function for 
one conjugate pair of bosonic valence quarks. The bosonic factor does not depend on the dynamical quarks and 
therefore goes beyond what is implied by the Toda lattice structure. This arises as a consequence of a singularity in 
the partition functions that contain a pair of conjugate bosonic quarks with the same mass. This singularity persists 
even for finite size matrices. In fact, it is present for an ensemble of 2 x 2 random matrices. Also the Toda lattice 
structure is valid for finite size matrices, and one can easily verify that the correct spectral density is obtained for an 
ensemble of 2 x 2 matrices. 

As a new result, we have also obtained explicit expressions for the spectral density when the absolute value of the 
fermion determinant is present in the QCD partition function. In this case a direct check from lattice simulations is 
possible for all values of the scaled chemical potential, as has already been done in the fully quenched case. 

The calculation of the unquenched microscopic spectral density has been a highly nontrivial test of the Toda lattice 
approach. We expect that it will be equally successful for all remaining (3 = 2 cases. For example, the method 
could be applied to the calculation of GUE S'-matrix fluctuations and parametric correlations [fj as well as to other 
non-Hermitian random matrix ensembles |72j . The extension of this approach to Dyson index (3=1 and (3 = 4 is still 
an open problem. 

The spectral density can be obtained from the replica limit of a Painleve equation that has been derived from 
fermionic partition functions only. On the other hand, if a family of partition functions satisfies a Painleve equation, 
a Toda lattice equation can be derived from the Backlund transformations of the Painleve system. Therefore, all 
systems in the same universality class have the same microscopic limit for the fermionic partition functions, the 
bosonic partition functions and the microscopic spectral density. 

Another important difference between Hermitian and non-Hermitian Gaussian random matrix theories is that the 
probability density of non-Hermitian theories do not factorize into the joint eigenvalue density and the distribution 
of the eigenvectors. The joint eigenvalue density is only obtained after a nontrivial integration over the similarity 
transformations that diagonalize the non-Hermitian matrices. This results in the appearance of a nontrivial universal 
factor in the joint eigenvalue distribution, in our case a modified if-Bessel function. The appearance of this function 
is a direct consequence of the chiral symmetry of the problem. 

The microscopic spectral density derived in this paper is valid for small temperatures and all values of the dimen- 
sionless parameters fj?F 2 V and rriEV . Within this range of parameters, where the sign problem develops and becomes 
severe, the effect on the spectral density can be followed analytically. A first manifestation of the sign problem is 
that the spectral density in the complex eigenvalue plane is no longer real. This occurs already in a parameter range 
where lattice simulations might be feasible. For larger values of ji 2 F 2 V the analytical expressions for the eigenvalue 
density can serve as tests for numerical methods that apply to the nonperturbative regime of QCD at nonzero baryon 
chemical potential. In the limit where fi 2 F 2 V ^> 1 the spectral "density" shows oscillations with a period of the order 
of the level spacing and an amplitude that diverges exponentially with /i 2 F 2 V. This behavior, which is not present 
in the quenched case, should be responsible for the breaking of chiral symmetry. This issue will be addressed in a 
future publication. 
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Appendix A. SOME DETAILS ON TWO FLAVOR AND GENERAL PARTITION FUNCTIONS 

For completeness we collect here the formulas for the remaining two partition functions with two flavors. As was 
observed in |38l |40| . partition functions that do not mix quarks or conjugate quarks in either the numerator or 
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denominator can be written solely as determinants of polynomials and Cauchy transforms. We find 



(x\y*;v) 



7 N f =0 




i n 



(z 2 - x 2 ) 

\(zf-y* 2 ) 



N,=0 



-1 



hN-i(y*) h N (y*) 
p N -i{x) p N {x) 



(78) 



and 



?N b =2 



(2/1,2/2;//) 



7 N f =0 
-■N 



(2/12/2)' 



N i 



r N -ir N -2(yl ~ yl){yiy2) 1 



N i= 

hjsi-2{yi) hff-i(yi) 
^-2(2/2) ^-1(2/2) 



(79) 



All remaining combinations of bosons and fermions and their conjugates can be obtained from the given ones by 
complex conjugation. 

Next we wish to comment on the relation between these formulas and in particular on the additional pole terms 
which occur in the two kernels containing Cauchy transforms, equations (|40(l and (|45(l . For random matrix models 
with real eigenvalues the kernel of the polynomials satisfies the Christoffel-Darboux identity, 



K N (s,t) 



1 p N (s)pN-l(t) -pjyft)piV-l(s) 
TN—X S 2 — t 2 



(80) 



For orthogonal polynomials in the complex plane Q28f) this relation is generally not satisfied so that (|35|l and l|38fl are 
different. In the Hermitian limit fi — > they are equal, due to the relation l|80|l . 

The Eqs. H39|l and (|79Jl or 144|l and 1|78J) can be related if the following Christoffel-Darboux identities are valid: 



A N -2(s,t) = 

Mv-i(M) = 



1 h N -i(s)tiN-2(t) - h N -i(t)h N - 2 {s) 



rjv-2 s 2 - t 2 

1 h N (s)p N -i(t) - p N (t)fiN-i(s) 



s,t £ R 



t 2 



s,teR 



(81) 



This is the case for Hermitian random matrix models, but in general there are no such Christoffel-Darboux identities 
for non-Hermitian random matrix theories. This identity follows from the observation that in general the orthogonal 
polynomials and their Cauchy transforms obey the same three-step recursion relation. Only the recursion involving 
the lowest Cauchy transform hk=o(y) is different from that involving the polynomial Pk=o{y), which results in the 
extra integral Q(s,t) over poles in l|4(J|) . or the single pole l/(s 2 — i 2 ) in l|45|l . respectively. Thus in the Hermitian 
limit l|39f) and (|79|l . and (|44|l and l|78|l become equal due to their Christoffel-Darboux identities 1)81(1. This explains 
the presence of pole terms in (|39H and l|44l) . In l|45|l this term can also be understood differently, as it ensures the 
correct normalization of the partition functions l|44|) to unity at equal arguments. 

Let us now turn to more general partition functions. We will briefly explain here how our generating function l|48|) 
follows from the results of [3!j (generalized to the chiral ensemble). There the following expectation value of complex 
eigenvalues is computed as 
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n 
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(82) 

In order to simplify the denominator the determinant can be expanded with respect to the next to last row. The 
remaining determinant only contains polynomials pk (mi ) and the kernel MN f -2 , y) ■ The invariance of determinants 
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under adding and subtracting rows and columns can be used to first eliminate all sums in the kernels A/jVf-2(^i! Tj)- 
Then all polynomials can be reduced to monomials, leading to the following Vandcrmonde like determinant 
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ml 



2(N f -2) 



r 5 9 



m 



N, 



2(N f - 



2(N f -l) 



(-1) 



N f 



(83) 



Here we have again made use of invariance properties of the determinant before arriving at the desired result. Sub- 
stituting this in (|H2J) we obtain (gBJ. 



Appendix B. THE MICROSCOPIC EIGENVALUE DENSITY WITH CONJUGATE QUARKS 

In this appendix we derive the microscopic eigenvalue density of the Dirac operator in a theory with conjugate 
fermionic quarks. Even if some or all of the fcrmionic flavors are conjugate quarks the Toda lattice equation holds. 
In the microscopic limit the partition functions 



r 

Z^' l ' n ({m f },{m c },{m*},z,z*;ti) = / [cL4]j, Y\_ det(D(/i) + m/) 

f =i 

i 

x Yl I det(D(n) + m c )\ 2 \ det(D(/i) + z)\ 2n e~ s ™, (84) 



satisfy the Toda lattice equation 

5 Z 8 Z * log Z u f " n ({mf}, {m c }, {to*}, z, z*; n) 

7771 (zz*) 2 Z » f ' l ' n+1 ({ m f}> { TO c}' i™*}: z i z *;t J -) z v f ' l ' n ~ 1 ({ m f}, {™c}, {m* c }, z, z*; fi) 
2 [Z^' l - n ({m f },{m c },{m* c },z,z*;v)] 2 

This follows by a direct extension of the argument given in section 2.3.2 of [2(j. In order to obtain the density with I 
pairs of conjugate quarks we consider the replica limit of the Toda lattice equation (|85|l with Nf = 

= f^HWa<),y^-y.(<>l^;.). (86) 

2 [Zl({m c },{m*};n)¥ 
Using H66f) we find the spectral density for a theory with / pairs of conjugate quarks 

(87) 

This density is positive and real for all masses and chemical potentials. This is fully consistent with having a real and 
positive measure in |0J. We have explicitly checked that this eigenvalue density is consistent with what we obtain 
using the orthogonal polynomial method as described in 

Finally, we give the most general result we have obtained using this method, namely the spectral density for a 
theory with Nf + I quarks and I conjugate quarks, 

p v f ' (z, z*, {to/}, {m c }, {m*}; fx) 

zz *Jli-,2 _2nTTu2 _2,2 7 „*. tA z u f ' l ' n ~ 1 ({mf},{m c },{m*},z,z*;fx) 
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